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Mechanical Assist Devices 
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• Motivation 


Over 5 million Americans and 20 million people worldwide suffer 
from Congestive Heart Failure (CHF) 

CHF patients are still treated with drug therapy, 
however, at late stage heart transplantation is traditionally 
the only treatment hope 




Mechanical Assist Devices 
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i Motivation 

- Need for assist devices is very high 

Need ! 26,000-60,000 / YR 

Donor hearts available ! 2,000-2,600 / YR 

(e.g. more than 4,000 patients were on the waiting list in 1999) 

- Need to find right match 

- Heart pump or VAD, for ventricular assist device, is being used 
as a temporary support to sick ventricle 
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Mechanical Assist Devices 
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• Motivation 

- VAD vs Drug treatment, 

recent study suggests that 

• Survival rate for VAD patients vs for patients receiving drug treatment 
After 1 year 52 % VS 24.7% (It al»o depend* on the method* and drug* iwed) 
After 2 years 22.9 % vs 8.1 % 

* Some patients who stayed in ICU because of short of breath can walk a 
block after 1 year assisted by VAD 



Mechanical Assist Devices 
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1 Motivation 

- VAD vs Drug treatment 

• However, complication rate for VAD is 2.35 times higher than that 
for drugs 

Complications include infections, bleeding, and mechanical malfunctions 
like motor failure, deformed tube and worn bearings 

Design improvements are needed to lower the risk, 

and possibly to use it as a permanent therapy (long-term device) 
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Mechanical Heart-Assist Devices 
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• Heart Valves 


• Ventricular Assist Device (VAD) 


- Piston Driven : Low speed. Bulky 

- Pneumatically Driven : Need external support equipment 

Rotary Pump 

- Axial Flow Pump : High speed. Small 
=> DeBakey VAD is based on this concept 


• Total Artif icial Heart 



Ventricular Assist Device 
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• Requirements 

- Simplicity and Reliability 

- Small size for ease of implantation 

- Supply 5 liter/min of blood against 100 mmHg pressure 

- High pumping efficiency to minimize power requirements 

- Minimum Hemolysis and Thrombus Formation 



Computational Issues 
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• Geometry / grid def inition 

Moving boundary 

• Solver 

Time accurate solver 

• Physical modeling 

Newtonian vs non-Newtonian 
Turbulence 

• Experimental & clinical data 


Solver: Viscous Incompressible Flow 
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• Formulation 

Can be viewed as a limiting case of compressible flow where the flow 
speed is insignificant compared to the speed of sound (Preconditioned 

compressible N-5 eq.) 

=> Artificial compressibility approach 
- Artificial Compressibility Method (Chorin, 1967) 

IN53D family of codes 
Merkle et. a! 
many more 

Or truly incompressible 

=o Pressure projection approach 

- MAC (Marlow and Welch, 1965) 

- Fractional Step Method (Chorin, 1968; Yanencko, 1971; Mar chuck, 1975....) 

- SIMPLE type Pressure Iteration (Caretto et a!., 1972; Patanka Si Spalding, 1972...) 
^ Use derived variables 

- Vorticity-Velocity (Fasel, 1976; Dennis et at, 1979; Hafez et at, 1988 

- Stream f unction-vorticity 



Artif icial Compressibility Method 
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• Formulation 

l^ + ^i=0 

- Introduces hyperbolic behavior into pressure field 

Speed of pressure wave depends on the artificial compressibility parameter, p. 

- The equations are to be marched in a time like fashion until the 
divergence of velocity converges to zero. 

=> Relaxes incompressibility requirement. 

Time variable during this process does not represent physical time step. 

For time -accurate solutions 

- Iterate the equations in pseudo-time level for each time step until 
incompressibility condition is satisfied 

=> Efficient sub-iteration is the key for success 


Artificial Compressibility Method (INS3D-UP) 
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> Time accuracy is achieved by sub iteration 

- Discretize the time term in momentum equations using second-order 
three-point backward-difference formula 

3 q"' - 4 q m + q"~ t (rhsY** 

- Introduce a pse Mrtime levcFand artificial compressibility, 

- Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 
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> Code I&fformance 

- Computing time : 60-120 ms/grid point /iteration (on £90 single epu) 

- Memory usage: Li no -relaxation 45 words/ grid point 

GMRES-ILU(0) 220 werds/grid point 



A pressure Projection AActhod(INS 3 D-FS) 


• Approach in generalized coordinates 

- Finite volume discretization 

- Accurate treatment of geometric quantities 

- Dependent variables - pressure and volume fluxes 

- Implicit time integration 

- Fractional step procedure 

Solve auxiliary velocity field first, 

then enforce incompressibility condition by solving a Poisson 
equation for pressure. 


e Code performance 

- Computing time : 80 ms/grid point/iteration (on C90 single epu) 

- Memory usage; 70 words/grid point 
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Pressure Projection Method 



• Fractional-step 

- Solve for the auxiliary velocity field, using implicit predictor step: 

-Jr («;-«;>= -vp m +/»(«,') 

At 

- The velocity field at time level (n+1) is obtained by using a 
correction step: 

T-(«r - *0= -Y 'p~ +/<0- V +Ku.) 

A t 

- The incompressibility condition is enforced by using a Poisson equation 

for pressure ( ) 

p'mp^-p 

vy-7-v- u 

* At 








A Current Challenges 
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• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important such as 

- Scalability (algorithmic & implementation) 

- Portability, transparent coding etc. 

• Computing resources 

- "Srid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 
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Parallel Implementation of INS3D 



MPI coarse grain + OpenMP fine grain 
24 zones / 2,8 Million points 


TEST CASE : SSME Impeller 
60 zones / 19.2 Million points 
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Parallel Implementation of INS3D 



Multi-Level Parallelism (MLP) 

INS3D-MLP : MLP routines + QpanMP 
Shared Memory MLP Organization for Origin 2000 


MLP Process 1 MLP Process 2 

Common /tocaV aa,bb Common nocai! aa,bb 



Common /global/ x,y ,2 




Tim® (tec) par Iteration 


Parallel Implementation of INS3D 
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INS3D-MLP (NAS MLP no pin-to-node) 
/ OpenMP 


TEST CASE : SSME Impeller 
60 zones / 19.2 Million points 



Number ol CPU* 


Number of CPU* 
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j/j Scripting for Acceleration of Grid Generation 
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INLET GUIDE VANES AND DIFFUSER 

Old IGV New IGV Old DIFF New DIFF 
No. of points (million) 7.1 1.1 8.0 1.6 

Time to build 1/2 day 10 see. 1/2 day 8 sec. 

Script timings on new grids based on SGI R12K 300 MHz processor 
Time to build script = 1 day for IGV, 1 day for DIFF 
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Parallel Implementation of INS3D 



RLV 2 nd Gen Turbo pump 

INS3D-MLP / 40 Groups 114 Zones / 34,3 M grid points 



Per processor Mflop is 
between 00-70. 

Cods optimization for 
oaoho based platforms 
Is ourrontly underway. 
Target Mflops Is to 
reaoh 120 per 
prooessor. 

Inoreaslng number of 
OpenMP threads Is also 
the main objeotlve for 
this effort. 




A SSAAB-rigl / Initial start flffl 

/tom fount to* 

After 2 1/2 rotation; 

• Status 

- 34.3 Million Points 

- 400 physical time steps in one rotation. 


- 

12 minutes wall time with 128 £PL/s on 
Origin platforms. One complete rotation 
requires 3.5-days wall-clock time with 
128 processors dedicated to the task. 

- I/O and memory management are critical 
for wall-clock time reduction 



• Issues / Needs 

- In reality, more than 10% of the 

super computing facility to one task is not 
always practical. 

- Need lOOx bigger supernode or use lower- 
fidelity method 

- Communication to/from designers and 
experimental group is a part of critical 

Particle Traces and Pressure 

technologies (in ©rid computing) 
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Example of Pulsatile Pump 



• Penn-State Artificial Heart 

Chimera ©rid for moving components 



This and other raeuft* wmr* first raportad by Kiri* at. al in 1991 ; 

•Computation of Incompressible Viscous Plows through Artificial Heart Daviess with Moving boundaries/ 
Proc. American Mathematical Society Summer Research on BlofMd Dynamics Conference, 

July 6-12, 1991, Seattle. W A 
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Example of Pulsatile Pump 



• Penn-5+a+e Artif icial Heart 

Analysis of time dependant data was an issue 








Issues in Axial flow VAD 



• Problems Related to Fluid Dynamics 
- Small size requires high rotational speed 
Highly eff icient pump design required 

-High shear regions in the pump may cause excessive bbod 
cell damage 

Minimize high shear regions 

-Local regions of recirculation may cause bbod cbtting 
Good wall washing necessary 

=>5mall size and delicate operating conditions make it difficult 
to quantify the flow characteristics experimentally 



A 

'AmsRe 


AmsRestJfrtCtster 


DeBakey VAD Development Timeline 



• Baseline Design 

1984 - NASA Johnson Space Center's David Saucier begins initial design 
work on axial pump VAD with Dr. DeBakey 

1988 - NASA/JSC and Baylor College of Medicine signs Memorandum of 
Understanding to develop the DeBakey VAD 

1992 - NASA/JSC begins funding the project 
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NASA/DeBokey VAD (Baseline Design) 



NASA /DeBakey Axial Flow VAD Impeller 



Rotational Speed : 12 ,600 RPM 
Flow Rate : 5 lit/min 




Traces Colored by Axial Velocity Magnitude 

"^OMO 0.285 0.610 


Rotational Speed :12,600RPM Flow Rate : 5 UtAnin 
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DeBakey VAD Development Timeline 


• CFD Assisted Design 

1993 - NASA/ARC is asked to develop CFD procedure to improve design and 
performance, D. Kwak and C. Kiris visit JSC to study the device 
The technology developed for rocket engine such as the Space Shuttle 
main engine was to be extended to blood f low simulation 



1994 . Kiris and K«ak begin work on 
design analysis using NAS 
supercomputers 

=> new design was proposed to 

INCLUDE AN INDUCER BETWEEN 
THE FLOW STRAIGHTNER AND 
THE IMPELLER 




Particle Traces Colored by Velocity Magnitude 






DeBakey VAD Development Timeline 


• CFD Assisted Design 



1994 -Kiris and Kwak continued design changes 

=» IMPROVE BEARING, NUB AND HUB EXTENSION DESIGN TO REDUCE BLOOD O.OTTING 




DeBakey VAD Development Timeline 



• Animal Tests 

1995 - Animal implantation* passed two-week requirements 

1996 - Full I design rights are granted to MicroMed, Inc. to produce the pump 

Began using bio -compatible titanium repbeing polycarbonate 

1997 - Configuration design finalized 


















DeBokey VAD Development Timeline 
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• Human Implantation in Europe 

1998 - On November 13, 1998, the first six DeBakey VADs are implanted in 
European patients by Roland Hetzer and DeBakey at the German 
Heart Institute of Berlin. One of the patients, fifty six year old Josef 
Pristov, Is able to return home and spend Christmas with his wife after 
a month's stay for recovery and monitoring at the clinic 
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NASA/DcBakey VAD 
Accomplishments to date (7/1/01) 



120+ patients implanted 

Number of patients currently ongoing with device 


Approved for 20 patients (14 male, 6 female) in a multi-center trial 
European trial 

Received "CE mark* (the EU equivalent to FDA approval) 

Results to date 

Fawrable compared to existing VADs 

Small incidence of thrombus is being investigated 
^ Further computational support is essential 



Summary and Discussion-1 



* Computational approach provides 

- a possibility of quantifying the flow characteristics: especially valuable 
for analyzing compact design with highly sensitive operating conditions 

- a tool for conceptual design and for design optimization 

* CFD + rocket engine technology has been applied 

- to modify the design of NASA/DeBakey VAD which enabled human 
implantation 

* Computing requirement is still large 

- Unsteady analysis of the entire system from natural heart to aorta 
Involves several hundred revolutions of the impeller 

- During on heart beat, impeller has 125 revolutions 

- With 1024 processors of Origin, one simulation (with several heart beat) 
from heart to aorta con be completed in months 





Summary and Discussion-2 
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* Further study is needed 

- to assess long term impact of mechanical VAD on human body, which 
requires modeling flexible wall and non-Newtonian effect and 
better downstream boundary conditions 

* There exist some geps between 

- CFD (assuming IT is a part of CFD applications) and 
biomedical expertise 



